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Topological configurations, monopoles and vortices, successfully describe quark confinement and 
the spontaneous breakdown of chiral symmetry. Despite their infinite action, these configurations 
are relevant due to a subtle cancellation between action and entropy. A natural explanation for 
this intrinsic fine-tuning is that smooth low action configurations exist which confine and which 
appear as singular topological objects in certain gauges. To reveal these confining semi-classical 
configurations, a new cooling method is proposed which largely reduces the action while preserv- 
ing the asymptotic quark-antiquark potential. First numerical results for a SU(2) gauge theory 
show that confining configurations with an average plaquette as high as 0.95 do exist. 
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1. Introduction 

Two low energy phenomena of Quantumchromodynamics (QCD) predominantly dictate the 
properties of matter under normal conditions: quark confinement and the spontaneous breakdown 
of chiral symmetry. It was already realised in the mid seventies that Yang-Mills theories admit 
topologically non-trivial configurations with finite action called instantons [1, 2, 3, 4]. It soon 
became clear that a dilute gas [2] or a liquid [4, 5] of instantons may explain the spontaneous 
breakdown of chiral symmetry. After the advent of computer studies of Yang-Mills theories by 
means of lattice gauge theories, configurations of low action have been constructed by cooling. 
It has been reported in that straightforward cooling of lattice configurations annihilates the string 
tension and therefore quark confinement [6]. Direct evidence of instanton generated confinement 
has not yet been established. 

The wide spread belief is that a different type of configurations which populate the Yang-Mills 
vacuum is responsible for quark confinement. For example, colour magnetic monopoles generically 
exist in certain gauges and might generate quark confinement be means of the dual Meissner ef- 
fect [7,8]. Alternatively, a vacuum populated with percolating vortices which carry flux in the cen- 
tre of the non-Abelian gauge group would equally well explain quark confinement[9, 10]. Equipped 
with rapidly increasing computational resources, these ideas were put to the test using lattice gauge 
simulations. It turned out that the early definitions of the centre vortices failed the quantitative 
test: the properties of those vortices were strongly dependent on the UV cutoff of the Yang-Mills 
theory. An important novel definition of the vortex matter was firstly proposed in [11, 12, 13], 
and the connection of those vortices has been established using lattice gauge simulations. It was 
pointed out for the first time in [14, 15] that these vortices are sensible in the continuum limit when 
the UV regulator of the theory is removed. In particular for the gauge group SU(2), this vortex 
texture of the vacuum explains many of the low energy features of Yang-Mills theory. In particular, 
deconfinement at finite temperatures appears as a vortex depercolation transition [16, 17, 18], and a 
close connection between quark confinement and spontaneous chiral symmetry breaking has been 
found: Removing vortices restores chiral symmetry as realise from the quark condensate [19], from 
the spectral density of the quark operator [20] or from an inspection of the Landau gauge quai^k 
propagator [21]. If the discretised version of the quark operator possesses acceptable chiral proper- 
ties even for singular vortex background fields, vortex only configurations also show spontaneous 
chiral symmetry breaking [22]. Up to now, the vortex picture for the gauge group SU(3) is less 
compelling: Centre projection falls short to reproduce the full string tension ('2/3 '-problem) [23], 
and the tight connection between confinement and spontaneous chiral symmetry breaking seems 
to be lifted [24]. The '2/3 '-problem is absent if the Laplacian Centre gauge is used to define the 
vortex matter [25]. As pointed out in [26, 23], however, the Laplacian vortex matter is not indepen- 
dent from the UV regulator. The practical consequences are that the vortex removed model theory 
has lost the property of asymptotic freedom and generically suffers from a low signal to noise ratio. 
Further SU(3) studies are needed to resolve these issues. 

Both, in SU(2) and SU(3), the scaling of vortex matter arising from the maximal centre gauge 
is a remarkable feature of the theory: the infinite action of the vortex world sheets is balanced by 
an equally large entropy [27, 28]. Since these vortices are defined in a particular gauge, i.e., the 
so-called maximal centre gauge, a natural explanation of this intrinsic fine tuning would be that 
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there is a smooth low action analogue of the vortex configurations which equally well induce con- 
finement. Unravelling these semi-classical configurations would facilitate an analytic first principle 
understanding of confinement. 

In this paper, we report a novel cooling technique which strongly reduces the action of the 
lattice configurations, but which preserves the long distance behaviour of the static quark anti- 
quark potential and thus quark confinement. 



2. Confinement preserving cooling 

Throughout this paper, we adopt the lattice regularisation of SU(2) Yang-Mills theory. Degrees 
of freedom are the matrices {x) associated with the link of the space-time lattice with lattice 
spacing a. We use the Wilson action 

•^wii = -i8 £ Ux\u^,{x)U,{x + pL)Ul{x + v)Ul{x)\. (2.1) 

For our cooling method below, we need to address the static quark-antiquark potential V (r) in terms 
of un-smeared Wilson loops. The expectation value of a rectangular Wilson loop W [r, t] with extent 
r in spatial and f » 1 /niQ in time direction is related to the static potential by: 



W[r,t] 



(2|a2)|'(r) exp{-V(r)f} , (2.2) 



where tuq is the mass gap of the theory, |2) is the state of two heavy, axial-dressed quarks [29, 30, 
31] and 1^2) is the true ground state in the quark-antiquark channel. For axial-dressed quarks, the 
overlap |(2|n2)p(r) is a small number which exponentially decreases with increasing r [30]. It is 
therefore strongly advisable to use Coulomb dressing or dressing via APE smearing, which both 
improve on the overlap, if the aim is the calculation of the static potential. 

The idea central to our novel cooling is to minimise the action (2.1), but to keep fixed the 
static potential V{r > r^). For our method, it will turn out below that fixing the potential for r = 
is sufficient to keep V{r> r^) unchanged. Hence, is a new scale which comes into play, and 
which is reminiscent of the so-called cooling radius employed in [32] for elaborate cooling. A 
straightforward but wrong idea would be to constrain the Wilson loop expectation value {W{r = 
rc,t)) during cooling. Since the size of the right hand side of (2.2) is dominated by the overlap 
rather than the exponential, this would merely imply that we freeze the poor overlap. This is not 
what we want. By contrast, cooling the lattice configuration should increase the overlap of the trial 
state with the true ground state. To this aim, we consider the ratio 



{W[r,t+l]) _ 
{W[r,t]) -"'^P 



{-V{r)a} . (2.3) 



A valid method would be to constrain the above ratio for r = rc. Since the numerator and the 
denominator consist of an average over many lattice configurations, this method is hardly feasible. 
To circumvent this practical problem, we firstly note that the Wilson loop W [r, t] [x] for a given 
lattice configuration depends on the position x on the lattice. The configuration average, however. 
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is position independent due to translation invariance ' . Consider now two Wilson loops at x and y 
with a separation much bigger than the inverse mass gap. We then find 

'W[r,t]{x) W[r,t]{y)) « (w[r,t]{x)) [w[r,t]{y)), jx-j] » -L . 

For pure Yang-Mills theory, the mass gap nic is set by the mass of the lightest glueball 
implying that t, := l/niQ is short ranged. This implies that Wilson loops separated by distances 
larger than are statistically independent. For large lattices, i.e., Na » we might therefore 
approximate the average over lattice configurations by an average over space time (and orientation; 
not shown) 

nrAi^)) « :^,"Lnr,m- (2.4) 

The total action which we are going to minimise is given by 

where A is a Lagrange multiplicator which ensures that 

L.W[r„f + l](x) 

during cooling. Any sufficiently large value for t should work. We have chosen 

t = rc 

throughout this paper. The action (2.5) can be minimised by standard means such as steepest decent 
method. The outline of the cooling approach is as follows: 

1 Generate a thermalised lattice configuration {Up}; 

2 Calculate the ratio R in (2.6); 

3 Calculate a change 8U of the link configuration which reduces the total action (2.5); 

4 Adjust the Lagrange multiplier A to ensure that r has not changed when you implement the 
change 8U ; 

5 Repeat steps 2-5 until the action (2.5) does not change anymore; 

6 Start again with step 1 to generate an new cooled configuration. 

In order to show that this approach is feasible, we here just report the change of the action 55 when 
the particular link Up{x) is changed to 5U while all other links are kept fixed: 

Bp{x)= Y,Uy{x + pi)Ul{x + v) Ul {x) . (2.7) 

Thereby, the (x) ai^e the well known staples, and [r^-, f] {x) is the sum of cut-open Wilson loops 
W[rc, f] starting at point x + pt and ending at point x. 

'in practical calculations, we sum over all positions to increase the signal-to-noise ratio. 
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Figure 1: Static quark-antiquark potential obtained from cooled configurations. 



3. Numerical results: 

First numerical results were obtained for a 16"* lattice for a gauge group SU(2) with j8 =2.4. 
The cooling radius was fixed to r<; = 3. We find that after 150 cooling steps further changes of the 
action are tiny. 

We firstly checked whether we recover the correct asymptotic behaviour of the static quark- 
antiquark potential from cooled configurations. To this aim, we evaluated the expectation value of 
rectangular Wilson loops averaged over cooled configuration. The would-be static potential V{r) 
is then obtained from 

-\n(w[r,t]) = V{r)t - function(r) , f > l/mc- 

\ / cooled 

For large values t > tf, this function is fitted to a linear function. If the ;t;^/dof is not satisfactory, a 
possible cure would be to increase tf. We studied the cases where the quark and the antiquark are 
aligned along the main cubic axis ((lOO)-direction) as well as along the diagonals (1 10) and (111). 
It turns out that the linear fit was never satisfactory for the orientation (100) of quark and antiquark 
while this problem did not occur for the diagonal directions. Further studies will be needed to 
trace out the origin of this effect. For the moment, we keep in mind that the cooled configurations 
might be plagued by sizable rotational symmetry breaking effects. Because of the constraint during 
cooling, we ensure that V{rc) = V(rc). The crucial question is whether V = V also holds for 
r > r^. Our numerical result obtained from 640 cooled configurations is shown in figure 1. The 
fitted string tension is given by 0.069(1) which compares quite well with value 0.075(2) from full 
configurations. After a sufficiently large amount of cooling steps (> 150 for my setup), the action 
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Figure 2: Action density obtained from cooled configurations with cooling radius = 3a, j3 = 2.4, 20^- 
lattice. 

density seems to stabilise. The constraint during cooling prevents the configurations to become 
the perturbative vacuum with trivial plaquettes only. An example of such a cooled configuration 
(cooling radius = 3a, p = 2.4) is shown in figure 2. We clearly observe a lump type structure in 
the vacuum. The overall achieved average plaquette of this configuration has been 0.95. 

4. Conclusions 

Over the last two decades, topological configurations, such as monopoles and vortices, have 
been successfully linked to the confinement phenomenon in Yang-Mills theories. It has been ob- 
served in the recent past that those degrees of freedom also describe the spontaneous breakdown 
of chiral symmetry. These configurations are relevant due to a subtle cancellation between action 
and entropy. A natural explanation for this intrinsic fine tuning is that this singular configurations 
have a low action counter part in a different gauge. In order to unravel these semi-classical con- 
figurations, a new type of cooling was introduced which largely reduces the action but preserves 
the asymptotic behaviour of the quark potential. The viability of the new cooling method has been 
explored for the SU(2) gauge theory for lattices as big as 20^*. A lump type vacuum structure has 
been observed using the action density. 

Future studies will address the questions: Does the action density of the cooled configura- 
tions scale like (a the lattice spacing) implying a finite action density in the continuum limit? 
What is the topological charge of the observed lumps, and what is their relation to instantons and 
calorons [33]? Does the observed structure match with monopoles and vortices in the correspond- 
ing gauges? 
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